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A method and flow meter for determining the flow rates of the components of a 
multiphase fluid 

The present invention relates to a method and flow meter for determining the composition and 
flow rates of individual components of a multiphase fluid, as defined in the preambles of claims 
1 and 12, respectively. 

The problem of how to meter oil-water-gas mixtures has been of interest to the petroleum 
industry since the early 1980s. Since then considerable research has been conducted into the 
development of a three-phase flow meter suitable for use in an industrial environment. 

There are several techniques and known instruments for measuring multiphase flow, as will be 
further described below. Such instruments need to be reasonably accurate (typically ± 5 % of 
rate for each phase), non-intrusive, reliable, flow regime independent, and suitable for use over 
the full component fraction range. In spite of the large number of solutions that have been 
proposed in recent years, no commercially available three phase flow meter yet meets all these 
requirements. 

The output of an oil/gas reservoir may vary greatly, depending on the location and age of the 
well. In addition to the oil and gas components, water, sand and wax may also be present in the 
produced well stream. Since the location and output of a well may vary so widely, the systems 
that have been designed to collect and process this output also vary considerably. The initial 
target of the oil industry to develop a universal multiphase flow meter to replace the traditional 
separation/single phase metering solution currently used, the fiscal monitoring of a well's 
output, have yet to be realised. 

Multiphase flow meters are increasingly used for well testing and allocation measurement. 
In order to optimise the production and life of an oil/gas field, operators need to be able to 
regularly monitor the output of each well in the field. The conventional way of doing this is to 
use a test separator. Test separators are expensive, occupy valuable space on a production 
platform, and require a long time to monitor each well because of the stabilised flow conditions 
required. In addition, test separators are only moderately accurate (typically ± 5 to 10 % of each 
phase flow rate) and cannot be used for continuous well monitoring. A three-phase flow meter 
could be used in the first instance instead of a test separator and in the long term as a permanent 
installation on each well. Such an arrangement would save the loss in production normally 
associated with well testing. Such loss is estimated to be approximately 2% for a typical 
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offshore installation. Allocation metering is needed when a common pipeline is used to 
transport the output from a number of wells owned by different companies to a processing 
facility. This is currently achieved by passing the output of each well through a test separator 
before entering the common pipeline. However, in addition to the disadvantages of the test 
separator described above, dedicated test pipelines to each well are also required. A 
permanently installed three-phase flow meter would offer significant advantages for allocation 
metering. 

According to a group of major oil companies, the accuracy requirements for a multiphase meter 
within a gas volume fraction range of 0-99% and water cut range of 0-90%, is 5-10% relative 
error on the liquid and gas flow rate and water cut measurement error within 2% abs. More 
accurate measurements were required for production allocation applications. Commercial three- 
phase flow meters are now generally capable of measuring individual phase fraction's flow rate 
to an uncertainty of less than 10% over a reasonably wide range of flow rates and phase 
fractions. There are two areas of operation which need further investigation if flow rate 
uncertainty is to be reduced still further using current combinational measurement techniques: 
flow regime dependency and individual phase velocity measurement. 

The present invention aims at providing a method and means which significantly reduce this 
uncertainty, particularly the uncertainty related to flow regime dependency. 

Some examples of commercially available non-intmsive multiphase meters such as those known 
from NO 304333, NO 304332, US 5,103,181, WO 00/45133 (figure 5) and US 6,097,786, 
measure the cross sectional composition and velocity of the phases to obtain flow rates. In order 
to provide accurate measurements, a homogeneous mixture in the cross section of the pipe is 
required. Effects due to inhomogenity in the longitudinal direction of the pipe is normally 
minimised by fast sampling of the cross-sectional composition. Multiphase meters are normally 
not mounted in a horizontal position due to the presence of laminar flow, where water is in the 
bottom of the pipe and gas at the top, which would distort the measurement. Consequently, to 
achieve homogeneous mixture in the cross section of the pipe of a multiphase meter, it is 
common practice to install the multiphase meters in such a way that the flow is flowing in an 
upward or downward direction. Laminar flow may then be avoided. However, when a 
multiphase mixture containing gas and liquid(s) are flowing in a vertical direction, annular flow 
can occur. Annular flow means that most of the liquid is distributed as a ring along the walls of 
the pipe and most of the gas is concentrated in the middle of the pipe. Annular flow distorts the 



wo 2005/057142 



3 



PCT/NO2004/000379 



measurement in a similar manner as laminar flow in a horizontal installation. In horizontal 
pipes pure annular flow where all the gas is in the middle of the pipe would normally only occur 
at higher gas fractions. However, when the flow is flowing in vertical pipes, severe 
concentration of gas in the middle of the pipe has been experienced even at medium flow rates 
(a few m/s) and gas fractions as low as 10%. Even a concentration of the gas in the middle of 
the pipe at lower gas fractions would introduce severe measurement errors. In practice, the 
liquid is rarely completely free of gas. In the context of this patent application, we define the 
degree of annular flow as, one minus the amoimt of free gas at the pipe wall divided by the 
amount of free gas in the middle of the pipe, as shown in the equation below. 

GVF 
GVF^ 

Where 

DOAF = Degree of Annular Flow 

GVFw = The amount of free gas at the pipe wall 

GWc = The amount of free gas in the middle of the pipe 

In other words, for a degree of annular flow of 1 .0 (or 100 %), all the gas volume is located in 
the middle of the pipe and all the liquid volume is distributed as a ring along the wall. 
Furthermore, for a degree of annular flow of 0 (or 0%), the gas volume is evenly distributed 
throughout the entire cross section of the pipe. Furthermore, for a degree of annular flow of 0.5 
(or 50%), there is twice as much gas volume in the middle of the pipe compared to the gas 
volume at the pipe wall. The values GVFw and GVFc are in the present invention only used in 
connection with caUbration of the mathematical model describing the relationship between the 
measured data and the degree of annular flow and not a direct part of the equations for 
calculation of the flow rates. 

NO 304333, US 5,103,181, US 6,097,786 and US 5,135,684 uses a nuclear densitometer. When 
a nuclear densitometer is used to measure the density, it is not possible to obtain fidl coverage 
of the cross section of the pipe. Hence, in order to obtain accurate measurements, it relies on a 
homogeneous mixture in the cross section. Typical commercial available nuclear detectors for 
density measurement, based on the Caesium 662 keV peak, has a circular area with a radius of 
2" and lower. For dual energy systems (x-ray and y-ray ) as described in US 5,135,684 and US 
6,097,786, the area is normally even smaller due to the need for a composite window in the pipe 
in order to allow radiation from the low energy x-ray radiation to go through the pipe. The cover 
area in a 2" pipe with a typical commercially available y-ray densitometer is typically 70-80 % 
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of the total cross sectional area of the pipe. However, when used in a 6'' pipe, it is difficult to 
achieve more than 30% coverage of the cross section of the pipe. One way to increase the 
coverage is to place the density measurement inside a venturi passage as in US 5,135,684. 
However, placing the nuclear density measurement inside a venturi passage also increases the 
amount of annular flow in the measurement section. When the source and detector is placed in 
the middle of the pipe, a too low density will be measured at annular flow. The error in the 
measurement will increase as the area of the pipe is increased. One way to compensate for this 
effect is to place the densitometer off-centre. However, the measurement errors due to gas 
concentration in the middle of the pipe would still be significant. 

Yet another way to minimise the effect of annular flow is to use a mixing device. US Re. 36,597 
describes a method where a positive displacement meter is used to both measure the total flow 
rate and homogenise the multiphase mixture in advance of the composition measurement. 
Aimular flow is then minimised; however, the multiphase meter becomes highly intrusive and 
fragile since it depends on a mechanical restricting or rotating device located in the multiphase 
stream. The repeatability of the measurement over time would also be vulnerable to sand 
erosion. Another way to reduce the presence of annular flow is to use a mixer. US 5,135,684 
refer to a method where a hold up-tank is used to homogenise the multiphase flow. However, 
the structure is highly intrusive, thus creating a pressure drop and hence limiting the production 
capabiUties from the wells. The performance of the mixer would also be dependent on the flow 
rate and pattern such as length of gas and liquid slugs and could therefore limit the operational 
envelope of such a multiphase meter. Another method based on mixing of the multiphase flow 
is described in US 6,272,934. 

Yet another way to reduce the effect of annular flow is to perform the composition measurement 
at the cross section of an annular venturi is shown in WOOO/45133, figure 1. However, this 
method is also intrusive and the repeatability of the measurement over time would also be 
vulnerable to sand erosion. 

Also known are multiphase composition and flow meters based on microwaves. 
US 4,458,524 discloses a multiphase flow meter that measures the permittivity (dielectric 
constant), density, temperature and pressure. Such a device uses phase shift between two 
receiving antennas to determine the permittivity. 

Other techniques are further known being based on resonance frequency measurement. 
Examples of such techniques are disclosed in WO 03/034051 and US 6,466,035. Techniques 
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based on resonance frequency measurements are normally limited to multiphase conditions 
where the loss inside the pipe is small and would therefore normally not work for high water cut 
applications and saline water due to the high dielectric loss of the mixture. US 5,103,1 81 
describes a method based on measurement of constructive and destructive interference pattems 
in the pipe. 

Also known are multiphase flow meters based on determination of permittivity of the 
multiphase mixture using capacitance and inductance sensors. Examples of such methods can be 
found in WO 00/45 1 33 and NO 304333, These devises utilise a lower frequency for 
measurements of permittivity compared to RF and microwave based systems, and are therefore 
much more sensitive to changes in the water salinity and droplet size of the multiphase flow. By 
using a lower frequency it is also difficult to design a single measurement unit capable of 
performing cross sectional measurements at both oil and water continuous flow conditions and 
any droplet size since the relative impedance of oil and water continuous flow differs by many 
orders of magnitude at low frequency. Capacitance and inductance based flowmeters are also 
more prone to measurement drift problems since the sensor and electrodes are a part of an 
electronic circuit. The parameters of the electronic circuit are measured by comparing them to 
reference values which need to be stable within a few picofarads in order to obtain the required 
accuracy. Such stability is difficult to achieve because of signal line capacitance drift, 
temperature drifts and stray capacitances in the system related to such as buildup of solids or 
oil/water film on the pipe walls. 

However, none of the above mentioned techniques are both non-intrusive and capable of 
performing accurate measurements of the flow rates at annular flow conditions or when the gas 
concentration is higher in the middle of the pipe. Thus, the main purpose of the present 
invention is to provide a method for accurate measurements of the flow rate of the individual 
components of a multiphase mixture at any flow regime, including annular flow and severe gas 
concentration in the middle of the pipe, without using a mechanical mixing or flow-conditioning 
device. 

It is a purpose of this invention to provide a method for identifying gas concentration in the 
middle of the pipe and to further compensate permittivity and density measurements for the 
measurement error related to the degree of annular flow. 

It is fixrther a purpose of this invention to provide an improved apparatus to avoid the above 
mentioned limitations of the performance of presently known techniques for multiphase flow 
measurements. 
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It is still further a purpose of this invention to provide a single low-cost structure for performing 
accurate measurement of oil, water and gas flow rates. 

And, it is a purpose of this invention to provide an almost non intrusive structure for performing 
the measurements without the need for an upstream mixing device. 

The method according to the present invention comprises the following steps: 

a. electromagnetic loss and phase measurements are performed in at least two 

directions of the pipe, 

b. the degree of annular flow is determined based on the measurements of step a, 

c. the permittivity of the flow mixture is calculated based on the results from 
steps a and b including correction for the degree of annular flow, 

d. the mixture density is measured and compensated for the degree of annular 
flow, 

e. the temperature and pressure are obtained, 

f. the velocity of liquid and gas are determined, and 

g. based on the knowledge of densities and permittivities of the components of 
the fluid mixture, and the result from the above steps a-f, the volume and mass 
flow rates of the gas and liquid or liquids of the fluid mixture are calculated. 

The flow meter according to the invention is further characterized by the features as defined in 
the independent claim 12. 

Dependent claims 2-11 and 13-21 define preferred embodiments of the invention. 

The invention will be further described in the following, by way of example, with reference to 
the figures, where: 

Fig. 1 shows a graphical representation of four examples of annular flow through a cross section 
of a 4" pipe. 

Fig. 2 shows a schematic longitudinal sectional view of a first flow meter according to the 
invention. 

Fig. 3 shows a schematic cross sectional view along the line IE - EI in figure 2, 
Fig. 4 is a sketch at a larger scale of a detail of figure 2, 

Fig. 5 shows phase difference as a function of frequency at low loss for the flow meter of figure 
2, 

Fig. 6 shows the electromagnetic field below the cut-off frequency TEn or at high loss inside a 
flow meter according to figure 2, 

Fig. 7 shows the electric field for the waveguide modes TEn and TMoi in the cross section of a 
flow meter according to figure 2, 
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Fig. 8 shows phase difference as a function of frequency at low loss for the flow meter of figure 
2, 

Fig. 9 shows the un-compensated frequency measurements for various liquids at pure annular 
flow for the flow meter of figure 2 

Fig. 10 shows measurements corresponding to figure 9, compensated for annular flow, 

Fig. 1 1 shows the un-compensated GVF (%Gas) measurements for various liquids at pure 
annular flow, 

Fig. 12 shows measurements corresponding to figure 11, compensated for annular flow. 
Fig. 13 shows a second embodiment of a flow meter according to the invention, 
Fig. 14 shows a third embodiment of a flow meter according to the invention, 
Fig. 15 shows a fourth embodiment of a flow meter according to the invention. 
Fig. 16 shows a fifth embodiment of a flow meter according to the invention, 
Fig. 17 shows a sixth embodiment of a flow meter according to the invention. 

The method according to the invention includes three main elements as follows: 

1) Measurements to determine the variation of gas concentration in the pipe cross- 
section, i.e. degree of annular flow. 

2) Measurements of the flow components, e.g. oil, water and gas, in a cross section of 
the pipe. This includes measurement of permittivity, density and obtaining the 
temperature and pressure of the flow. The permittivity and density measurements are 
compensated for the variation of gas concentration in the cross section of the pipe. 

3) Measurements of liquid and gas velocity. By combining the measurements from 
points 2) and 3) above and knowing the cross-sectional area of the pipe (sensor) and 
the density of oil, water and gas - it is possible to calculate the oil, water and gas 
volumetric and mass based flow rates 

The weakness of existing multiphase meters is mainly related to two factors, namely that: 

1) Existing multiphase meters rely on a homogenous mixture between oil, water and gas in 
the cross section of the pipe. At large variation of gas concentration in the pipe cross- 
section, great measurement errors will occur. In Fig. 1 is shown a graphical 
representation of the effect on the GVF (% Gas) measurement with a typical gamma 
densitometer based on a 2" detector 1 and a gamma source 34 for four examples of 
annular flow where all the gas 37 is in the middle of the pipe and all the liquid 38 is 
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distributed along the pipe wall. Although this is an extreme case since all the gas is 
concentrated in the middle of the pipe, it illustrates that large errors will occur in the 
measurements. 

2) Need of using a mechanical mixing element. Some multiphase meters use a mechanical 
mixing device to homogenise the multiphase flow. A mixer would reduce the amount of 
annular flow; however, it makes the meter highly intrusive. Some of the mixers may 
also contain moving mechanical objects that are vulnerable to sand erosion and could 
even be damaged by fast variations in the velocity associated with start up of a well. 

The imiqueness of the present invention is the ability to detect presence and degree of annular 
flow and compensate the measurement error related to the degree of aimular flow. 

The attached Fig. 2 illustrates a flow meter according to the invention. The antennas 16, 17, 18, 
19, 20 and 21 , as can be seen in more detail in Fig 4, are in effect coaxial conductors that are 
inserted into the pipe designed such that the centre conducting wire 22 which is isolated from 
the pipe wall 24 by a dielectric material 23 such as plastic or ceramic. Three of the antennas are 
in this example used as transmitters, and are therefore given prefix a Tx, and three of the 
antennas are used as receivers and are therefore given a prefix Rx. The direction of the flow is 
illustrated by an arrow 25. A densitometer 29 is also a part of the flow meter. An electronic 
system capable of transmitting and receiving a broadband signal (typical 10 Mhz - 4,0 Ghz) on 
the various antennas, a computer and devices for measurement of temperature and pressure are 
also parts of the flow meter. But for simplicity, these devices are omitted from all the figures 
and the fiirther discussion since it will be clear to the skilled person how they may be 
implemented. The pipe arrangement of figure 2 may also be referred to as a sensor in the 
context of this patent application. 

Loss and phase measurements are performed by measuring the received power and phase 
difference of a broad band signal (typical 10 Mhz - 4,0 Ghz) which is transmitted from a 
sending antenna and received at two receiving antennas located at a different length from the 
sending antenna. The measurement is done at at least two and preferable three planes in the pipe 
where one plane is at the cross section, a second is at the longitudinal direction and a third at an 
angle (such as 45 degrees) to the flow direction. The frequency is typically varied form 10 Mhz 
until 4,0 Ghz depending on the pipe diameter. By recording the frequency at at least three 
predetermined phase differences and using a calibration constant for the system, the permittivity 
within the pipe can be measured in all directions. Based on the measurements of the permittivity 
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in the various directions, the degree of annular flow can be measured and compensated for by 
using a mathematical model such as neural networks since the various measurements are 
differently affected by concentration of gas in the middle of the pipe. 

Figure 5 shows the phase vs. frequency measurement of a flow meter according to figure 2 at 
high loss inside the pipe. When the loss inside the pipe is high, the electromagnetic field will 
behave according to plane wave theory. When the flow is well mixed, the phase vs. fi^equency 
curve 4 is more or less linear as shown in figure 5. When gas is concentrated in the middle of 
the pipe, the curve 5 is more non-linear. Hence, by analysing the distribution of the phase 
measurements, annular flow can also be detected and compensated for. One way to develop a 
mathematical model describing the behavioxir of this feature is to use empirical data and train a 
neural network to detect annular flow and compensate the measurement. The network could be 
trained to predict the corresponding well mixed permittivity or measurement fi*equency, where 
well mixed permittivity or fi-equency in this context means the theoretical permittivity or 
firequency that would have been measured in an equivalent homogeneous multiphase mixture. 

Laboratory tests based on this method has provided 100% identification of annular flow, and 
large reduction in measurement errors due to annular flow. 

When the pipe acts as a wave guide (low loss), the permittivity is measured by measuring a 
fi*equency related to the cut-off fi-equency for some of the wave-guide modes of the pipe. The 
distribution of the electric and magnetic field within a pipe is firequency dependent. Table 8.9, 
on page 425 of Fields and Waves in Communication Electronics written by S. Ramo, J.R 
Whinnery and T. V. Duzer (John Wiely & Sons, 1964, second edition) shows the electric and 
magnetic field lines and the corresponding cut-off fi-equencies for the various TM and TE 
classes of waves (waveguide modes) of a circular pipe. Figure 6 shows the electromagnetic field 
inside a flow meter according to fig. 2 at high loss or at a fi-equency well below the cut-off 
fi-equency TEj i of the pipe. Below the lowest cut-off fi-equency of the pipe, which is TEn, the 
electromagnetic field will propagate according to plane wave theory as illustrated by the arrows 
6 in the Fig. 6. Fig 7 shows the electric field lines of the two lowest wave-guide modes of a 
circular wave-guide, namely TEi i 7 and TMoi 8. When the electric field (E-field) in the pipe 
changes firom plane wave propagation into TEn, a step occurs in the phase difference 9 of the 
receiving antennas as shown in figure 8. By applying a fi-equency sweep on the transmitter and 
measuring the firequency at at least three predetermined phase differences, the frequency of the 
phase step, which is a measure of the cut-off firequency TEu of the pipe, can be calculated. 
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A measurement of the permittivity within the pipe performed in the longitudinal direction of the 
pipe can be obtained by placing a microwave reflector 12, such as a cross or fin with a length of 
approximately 0.5 pipe diameters, at a predetermined distance from the transmitting antenna. 
Such an arrangement is shown in figure 13 where a microwave reflector 12 is placed upstream 
the transmitting antenna Tx3 19. 

Figure 8 shows the measured phase difference vs. frequency for the sensor arrangement of 
figure 15 at low loss. The frequency location of the phase change 9 corresponds to the cut-off 
frequency of the wave-guide mode TEn of the venturi throat 10. The frequency location of the 
phase change 13 of figure 8 corresponds to the first half-wave reflection between the 
transmitting antenna Tx3 19 and the reflector 12, which is above the cut-off frequency TEn of 
the large pipe 1 1 and is a function of the permittivity inside the pipe. Similarly, the phase 
change 15 corresponds to the frequency of the first half-wave reflection between the 
transmitting antenna Tx3 19 and the reflector 12 which is above the cut-off frequency TMoi of 
the large pipe 1 1 , The frequency locations of the phase shifts 9, 13 and 15 are also a function of 
the permittivity within the pipe. Hence, by measuring the frequency location of the phase shifts 
9, 13 and 15 and using a calibration constant for the system, three measurements of the 
permittivity within the pipe can be obtained. TEn 7 and TMoi 8 have different E-field 
distributions in the cross section of the pipe as shown in figure 7. Since one measurement is 
performed in the cross section and two measurements are performed in the longitudinally 
direction based on the waveguide modes TEn and TMoi, the three measurements are differently 
affected by an inhomogeneous mixture in the cross section of the pipe, such as gas 
concentration in the middle of the pipe (annular flow). The measurements can be obtained by 
performing a frequency sweep on Txl 16 and measuring the frequency at at least three 
predetermined phase differences on Rxl 17 and Rx2 18. The frequency location of these phase 
shifts are used as a starting point for a second and third frequency sweep transmitting on Tx3 19 
and recording the phase between Rx3 20 and Rx2 18. First, the frequency is further increased 
recording the frequency at at least three pre-determined phase differences in order to identify the 
frequency location of the phase shift 15. Then the frequency is reduced while recording the 
frequency at at least three predetermined phase differences for identifying the frequency 
location of the phase shift 13. Based on the recorded frequencies for the phase shifts 9, 13 and 
15 and using calibration constants for the system, the permittivity within the pipe can be 
measured. Based on these three measurements of the permittivity, the degree of annular flow 
can be calculated using a mathematical model such as neural networks since the three 
measiu-ements are differently affected by the degree of annular flow. The model can be derived 
based on empirical data obtained by measuring the phase shifts 9, 13 and 15 for a wide range 
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known degrees of annular flow. The calculated degree of annular flow is then used to 
compensate the permittivity and density measurement using an experimental derived 
mathematical model such as a neural network. 

At low loss, the degree of annular flow can also be measured by selecting a measurement 
frequency that is well below the cut-off frequency TEn of the pipe and measuring the loss in 
two or more planes of the pipe. When using this method at low loss to measure the degree of 
annular flow, the reflector 12 can be omitted from the flow meter, making it less intrusive. Since 
the cut-off frequency of TEn 9 is a function of the permittivity of the multiphase mixture inside 
the pipe, the measurement frequency will vary as a function of the permittivity. At the 
measurement frequency, the pipe would not act as a wave guide and consequently the E-field 
will be as shown in 7 of Fig. 7. At a frequency that is well below the cut off frequency of TEn, 
the E-field will propagate according to plane wave theory as shown in 6 of figure 6. The loss is 
measured by transmitting on the selected measurement frequency well below the cut-off 
frequency of TEn on antenna Txl 16 of figure 2 and recording the received power on antenna 
Rxl 17 and antenna Rx2 18 of figure 2. Then, using the same frequency, the next step is 
transmitting on antenna Tx3 1 9 and receiving on antenna Rx2 1 8 and Rx3 20. A third plane can 
be obtained by transmitting on antenna Tx3 19 and receiving on antenna Rx2 1 8 and Rxl 1 7. 
Based on the measured loss in the two or three measurement planes, the degree of annular flow 
can be calculated using an experimentally derived mathematical model and used to compensate 
the permittivity and density measurements using a mathematical model such as a neural 
network. 

Figure 9 shows the measured frequency vs. the theoretical frequency based on the cross 
sectional measurement plane for a gas void range (GVF) of 0-64% at a degree of annular flow 
of 1 .0 and a wide range of liquids. The measured frequency is the average frequency of three 
different predetermined phase differences and the theoretical frequency in this context means 
the theoretical value that would have been measured in an equivalent homogeneous multiphase 
mixture where the degree of annular flow is 0.0. Figure 10 shows the compensated frequency 
measurements based on the method described above using a neural network as the mathematical 
model for calculating the theoretical homogeneous frequency value. Figure 1 1 shows the 
uncompensated GVF (%Gas) measurements for a GVF range from 0-64% at a degree of annular 
flow of 1 .0 and a wide range of liquids. Figure 12 shows the measurements of the GVF 
compensated for the degree of annular flow. 
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In order to calciilate the oil, water and gas fractions in the cross section of the pipe, a 
measiarement of the cross sectional density is also required. This measurement would also be 
affected by the degree of annular flow. Knowing the degree of annular flow in the pipe, the 
density measurement can be compensated for the effect on the measurements in a similar 
manner using an experimentally derived mathematical model, such as a neural network. 

The density measurement is, according to the invention, performed in two ways depending on 
the application: 

D Gamma rav absorption (figure 13 and 17). B y measuring the gamma ray absorption 33 
of the multiphase mixture based on radiation from a gamma source 34 and knowing the 
absorption coefflcient of oil, water and gas and the permittivity of the multiphase 
mixture and the permittivity of oil, water and gas, it is possible to calculate the mixture 
density in an iterative calculation. As a part of this iteration, the gamma ray absorption 
measurement can be corrected for the presence of annular flow by a mathematical 
model, such as a neural network. 

2} Venturi mass flow measurement (figure 14 and ISV A venturi can be used to measure 
the density of the mixture. The pressure drop across the inlet of a venturi 30 is a 
function of the mass flow and density of the multiphase mixture. Furthermore, the 
pressure drop across the outlet of the venturi 31 is a function of the mass flow, density 
and compressibility of the multiphase mixture. Combining the pressure measurement of 
the inlet 30 and outlet 3 1 of the venturi, together with the measurement of the gas and 
liquid velocity from cross correlation (described in section below), it is possible to 
calculate the mixture density in an iterative fashion. However, at a degree of annular 
flow above 0, an error will be introduced to the density measxirement. Then, as a part of 
the iteration, the measurement can be corrected for the degree of of annular. 

A combination of a gamma ray absorption (pt 1) and venturi (pt 2) measurements, e.g. as shown 
in figure 16 and 17, may also be used. This combination can in some cases extend the operation 
envelope of the measurement system and enhance the measurement accuracy. In figurel6 and 
1 7 the gamma densitometer is placed inside the venturi throat together with the antennas such 
that the measurements can be performed under the same conditions. However, the gamma 
densitometer 33 and 34 may also be placed at either end of the sensor, but then it requires a 
compensation model to correct for the difference between the venturi throat 10 and pipe 1 1 . 
This correction model can be derived based on empirical data. When the antennas are placed 
inside a venturi throat, the antennas 16, 17 and 1 8 has to be approximately 0.5 throat 10 



wo 2005/057142 



13 



PCT/NO2004/000379 



diameters from the start of the venturi throat; otherwise, the measurement of the cut-off 
frequency will be influenced by the diameter of the large pipe 1 1 . Alternatively, in order to 
make the sensor more compact, the cross sectional measurement in such an arrangement can be 
placed more in the middle of the venturi throat, as shown in figure 17. Two additional antennas 
are now added to the venturi throat, namely Rx4 35 and Rx5 36. The cross sectional 
measurement can now be obtained by transmitting on Tx3 1 9 and receiving on Rx4 35 and Rx5 



The sensor is used to measure the composition and velocity (liquid and gas) of the multiphase 
mixture. Below is a more detailed description of the equations involved. 

Velocity Measurement 

By continuously transmitting and measuring the loss at the antenna pair Txl 16 / Rx2 1 8 and 
Tx221 / Rx3 20 located at a known distance S+L 26, one can create two time varying signals 
that are shifted in time equal to the time it takes the multiphase flow to travel between the two 
antenna pairs. The measurement frequency is selected such that little energy is going in the 
longitudinal direction. At low loss the frequency would typical be substantially below the cut- 
off frequency TEn of the pipe. By cross correlating the two signals using the formula: 

Equation 1: 



where x(t) and y(t) are the sampled signals, the time delay x can be calculated. The time delay x 
between the signals x(t) and y(t) is a measure of the time it takes a disturbance in the flow to go 
from the first to the second pair of antennas. Using high frequency signals to measure the flow 
disturbances also enables use of high sampling rates since it is possible to perform single 
measurements within a few micro seconds. Hence the signal contains information about small 
variations such as small gas bubbles in the liquid phase or water droplets in the oil phase or oil 
droplets in the water phase that typical represents the velocity of the liquid, and large variations 
such as gas slugs that represent the velocity of the gas phase. By applying the appropriate 
filtering of the sampled data and statistically sorting of the cross correlated velocities, it is 
possible to obtain a measure of both the liquid and gas velocity (viiq and Vgas). 
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Composition Measurement 

However, in order to measure the flowrates of oil, water and gas, it is required to measure the 
cross sectional composition (%oil, %water & %gas) of the multiphase mixture of oil, water and 
gas. By measuring the mixture permittivity Smix and mixture density Pmix the following equations 
can be used: 

Equation 2 : 

<I>oil+<I>water+<l>gas= 1 

where : 

^oii = Cross sectional volume fraction of oil 
^water= Cross scctioual volume fraction of water 
^gas = Cross sectional volxmie fraction of gas 

Equation 3: 

Pwater "'"^gas X pgas ~ Pmix 

where: 

Poii = Density of oil 
Pwater ~ Density of watcr 
Pgas = Density of gas 
Pmix = Measured density 

A temperature and presswe measurement is also required in order to compensate the above 
density parameters for temperature and pressure variations but, for simplicity, these will be 
omitted in the following discussions of the measurement principle. 

The Bruggeman mixing equation relates the permittivity (dielectric constant) of a two 
component mixture to the volxmie fractions of the components. If the two component mixture is 
droplets as an inner phase dispersed in a continuous medium of an outer phase, the equation 
become: 
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Equation 4: 



where: 



^inner ^mix 

F — £ 

^ inner outer 



1 

>i3 



dinner 
Souter 

<I>mner 



= 1-- 



mner ^ outer 



= Permittivity of the inner phase (dispersed phase) 
= Permittivity of the outer phase (continuous phase) 
= Measured permittivity of the mixture 
= Volume firaction of inner phase (dispersed phase) 
= Volume fraction of outer phase (continuous phase) 



A temperature and pressure measurement is also required in order to compensate the above 
permittivity parameters for temperature and pressure variations but, for simplicity, these will be 
ignored for the following discussions of the measurement principle. 

The equation above can also be used for a three-phase mixture such as oil, water and gas in 
which the inner phase is a well mixed combination of two of the phases dispersed in an outer 
phase. E.g., an inner oil/water mixture may be dispersed in an outer continuous media of gas 
and similarly, gas bubbles may be dispersed in an outer continuous media of an oil/water 
mixture. 



The lowest cut-off frequency of a circular wave guide, such as the pipe section of the flow 
meter, is TEn at: 



Equation 5: 



^ ^ 0.293 

where: 

fc = Cut-off frequency 

r = Radius of pipe 

e = Permittivity (dielectric constant) inside the wave guide 

(pipe) 

|Li = Permeability inside the wave guide (pipe) 
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Below the cut-off frequency, the electric field 6 will propagate according to plane wave theory 
as illustrated in figure 6. At low loss in the pipe and above the cut-off frequency f^, the electric 
field in the pipe is shown in 7 of fig 7 which correspond to TEn. When the field in the pipe 
changes from plane wave propagation into TEn, a step occurs in the ph£ise difference of the 
receiving antennas Rxl 17 and Rx2 18 of Fig. 2. The phase step is illustrated in 9 of Fig. 8. By 
applying a frequency sweep on the transmitter Txl 16 and measuring the frequency at at least 
three predetermined phase differences between the two receiving antennas, the frequency 
location (measured frequency) of the step change in the phase difference between the receiving 
antennas can be measured. Then, the measured frequency is a measure of the cut-off frequency 
fc of the pipe. 

Eqxxation 5 can be rearranged as: 
Equation 6: 




where: 

0.293 




fc = Frequency of electromagnetic wave (cut-off frequency 
ofTEu) 

8 = Permittivity (dielectric constant) inside the pipe 

hence k2 can be determined by measuring the frequency fc with a known permittivity inside the 
pipe such as vacuum where the permittivity is 1,0. 

The permittivity of the mixture at high loss inside the pipe (sensor) is measured by applying a 
frequency sweep to one of the transmitting antennas 16 or 19 and recording the frequency at at 
least three predetermined phase differences between two of the receiving antennas 17/18 or 18 
/ 20 located at a distance S and distance L from the transmitting antenna. Below the cut-off 
frequency or when the loss inside the pipe is large, the electric field will propagate according to 
plane wave theory. The phase difference between the two receiving antennas represents the 
wave travel time between the two points and can be written as: 
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Equation 7: 



where: 



AS = A 

AS 
A0 

antennas 
X 



A6 
In 



L-S (26) 

Phase difference between receiving 
Wavelength 



According to plane wave theory, the velocity of an electromagnetic wave can be expressed as: 



Equation 8: 



where: 



f 

£ 

c 



Frequency of electromagnetic wave 
Wavelength of electromagnetic wave 
Permittivity (dielectric constant) inside the pipe 
Permeability inside the pipe 
Speed of light 



Since the frequency is measured at predetermined phase difference, equation 6 and 7 can be 
combined giving: 



Equation 9: 



where: 



k, = 



cA6 



f 

8 



Frequency of electromagnetic wave 
Permittivity (dielectric constant) inside the pipe 
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ki can be determined by measuring the frequency at the phase difference A9 with a known 
permittivity inside the pipe. 

The permittivity within the pipe is measured in at least two directions. First, the transmitter is 
sending on Txl 16 and receiving on Rx2 1 8 and Rxl 1 7 (figure 2) performing a measurement of 
the permittivity in the cross section of the pipe. Then the transmitter is sending on Tx3 19 and 
measuring on Rx2 1 8 and Rx3 20 performing a measvirement of the permittivity in the 
longitudinal direction of the pipe. It is also possible to perform measurements of the permittivity 
by sending on Tx3 19 and receiving on Rxl 17 and Rx2 18 and hence performing a 
measurement that lies between the cross section and longitudinal measurements. 

The effect of annular flow on various measurement directions may be explained as follows. 
When the flow is well mixed the phase-difference vs. frequency would be almost linear. If the 
flow is annular, which distorts the symmetry of the L and S path from the transmitter to the 
receivers, the phase difference would be much more curved. The longitudinal antennas (18, 19 
and 20), as shown in Fig. 2, are less affected by annular flow since the symmetry is maintained 
also at annular flow. 

By measuring the frequency at several predetermined phase differences, it is possible to both 
detect and compensate for the effect on the measurements. Experimental data has shown that the 
effect on the measurement is related to the slope (d9/df) of the phase difference. One way to 
compensate for the error introduced by annular flow is to first train a nexiral network to calculate 
the degree of annular flow. Then a second neural network could be trained to compensate for 
the error in the permittivity measurement related to the degree of annular flow. 

The presence of annular flow can also be measured by measuring the loss in the longitudinal 
and cross-sectional direction. First, the transmitter is sending on Txl 16 and receiving on Rx2 
1 8 and Rxl 17, thus performing a measurement of the relative loss in the cross section of the 
pipe. Then the transmitter is sending on Tx3 1 9 and measuring on Rx2 1 8 and Rx3 20, 
performing a measurement of the relative loss in the longitudinal direction of the pipe. It is also 
possible to perform measurements by sending on Tx3 1 9 and receiving on Rxl 1 7 and Rx2 1 8 
and hence performing a measurement that lies between the cross section and longitudinal 
measurements. At a degree of annular flow above 0, the longitudinal measurement would be 
different compared to the cross-sectional measurement. The measurement has to be made in 
such a way that the pipe does not act as a wave guide. One way to achieve this is by selecting a 
measurement frequency that is below the measured cut-off frequency for TEu. 
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Yet another way to obtain measurement of the permittivity in the longitudinal direction of the 
pipe is to use a sensor arrangement as shown in figure 15 and described in the foregoing section. 
In this case the calibration constant k2 of equation 6 would be equal to the fi-equency location of 
the phase shift 13 and 15 of figure 8 with vacuimi inside the pipe (sensor). 



Measurement of gamma ray absorption is a widely used technique for density measurement. 
This technique takes into account that absorption of photon beam radiation in any material in 
the pipe (flow meter) can be expressed by the formula: 



Equation 10: 

where: 

No = Empty pipe count rate (radiation) 

N = Measured count rate (radiation) 

= Radiation mass absorption coefficient of the material inside the pipe, 

d = Transmission length of the radiation through the cross-section of the 

pipe 

p = Density of the material inside the pipe 



By measuring the count rate with a media inside the pipe with a known absorption coefficient 
such as fi-esh water, the parameter d can be determined according to equation 1 1 : 



Equation 11: 



N 

I ✓ fresh-water ^ 

N 



P fresh —water fresh -water 



where : 

No = Empty pipe count rate (radiation) 

Nfi«sh-water ~ Measured count rate (radiation) in fresh water 

M-fresh-water = Radiation mass absorption coefficient of firesh water 

pfi^sh-water = Dcusity of fi-csh-watcr 
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The density measurement does not cover the whole cross-sectional area of the pipe, 
consequently it relies on a homogeneous mixture in the cross section. The cover area in a 2" 
pipe with a typical commercial available y-ray detector is typical 70-80 % of the cross-section. 
However, when used in a 6" pipe, it is difficult to achieve more than 30% coverage of the pipe 
cross section. Nevertheless, knowing the degree of annular flow in the middle of the pipe, it is 
possible to compensate the measurement to provide a more correct measurement of the cross- 
sectional liquid and gas ratio. The compensation algorithm can either be derived from a 
geometrical description of the nuclear coverage area inside the pipe or by using an 
experimentally derived mathematical model, such as a neural network trained to correct the 
measurements. 

Yet another way to measure the density is to use a venturi mass flow meter as shown in figure 
14 and figure 15. Any restriction in the pipe will result in a change in the velocity of the 
multiphase mixture and introduce a pressure drop across the restriction. Based on the theory of 
fluid dynamics, the square root of the pressure drop 30 is proportional to the total mass flow rate 
in the pipe. A venturi tube is a structure where the pipe diameter is gradually reduced into a 
section of the pipe with a smaller diameter. The smaller section may be short or a relative long 
section. Then the diameter is gradually expanded to the original size of the pipe. Mass flow 
measurements with such a structure are described in ISO 5167-1 and ISO 5167-4. 

According to ISO 5 1 67-1 , the mass flow rate can be calculated as: 

Equation 12 : 



Qm = 



C 71 



where: 



Qm 



Total mass flow rate 



Ap 



C 



d 



P 



Discharge coefficient 

Diameter ratio between venturi throat and pipe 
Diameter of venturi throat 

Measured pressure drop between inlet and venturi throat 
Density of the multiphase mixture 
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The pressure recovery at the outlet of the venturi will mainly depend on the mass flow rate, 
density, compressibility and viscosity of the multiphase fluid and the length and roughness of 
the venturi throat 10. When the gas content of the multiphase mixture is high, the pressure 
recovery at the outlet of the venturi will be greater compared to a multiphase mixture with low 
gas content. Hence, by combining equation 12 with a measurement of the pressure recovery at 
the outlet of the venturi 3 1 , it is possible to obtain a measurement of the density of the 
multiphase mixture. 

Equations 1-12, together with the correction functions, are typically solved in an iterative 
fashion to derive the oil, water and gas volumetric and mass based flow rates, using a computer 
as an integrated part of the flow meter. 

Although several flow meters for measuring the flow rates of oil, water and gas have been 
described as examples for utilizing the present invention, the invention may also be used in 
other areas such as measurement of multiphase slurries with air or gas within the process 
industry, measurement of multiphase streams of coal and air or steam and water within the 
power generation industry. Furthermore, it will be clear to the skilled person that the invention 
is not limited to the embodiments described herein, but may be varied and modified within the 
firame of the invention defined by the features set forth in the appended claims and the 
equivalents thereof. 



